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Thermodynamic length is a path function that generalizes the notion of length to the surface 
of thermodynamic states. Here, we show how to measure thermodynamic length in far-from- 
equilibrium single molecule experiments using the work fluctuation relations. For these microscopic 
systems, it proves necessary to define the thermodynamic length in terms of the Fisher information. 
Consequently, the thermodynamic length can be directly related to the magnitude of fluctuations 
about equilibrium. The work fluctuation relations link the work and the free energy change during 
an external perturbation on a system. We use this result to determine equilibrium averages at inter- 
mediate points of the protocol in which the system is out-of-equilibrium. This allows us to extend 
Bennett's method to determine the potential of mean force, as well as the thermodynamic length, 
in single molecule experiments. 

PACS numbers: 05.70.Ln, 05.40.-a 



Modern experimental techniques allow the manipula- 
tion of single molecules and the measurement of the ther- 
modynamic properties of microscopic systems [l], 0, S 0, 
5]. For example, Collin et al. Q recently measured the 
work performed on a single RNA hairpin as it was folded 
and unfolded using optical tweezers. From these out- 
of-equilibrium measurements, they extracted the equi- 
librium free energy change using the recently discovered 
work fluctuation relations [1 0, H- [Eq. (JUJ] These re- 
lations, which connect the free energy change and the 
work done on a system by an external perturbation, re- 
main valid no matter how far the system is driven away 
from thermal equilibrium. 

In this Letter, we will demonstrate that free energy is 
not the only important quantity that can be extracted 
from out-of-equilibrium work measurements; we can also 
measure the thermodynamic length [tj [l^, 11, 12, [l3|, 
14, 15, [ill. Thermodynamic length is a path function 
that measures the distance along a path in thermody- 
namic state space. This is in contrast to the free energy 
change, a state function which depends only on the ini- 
tial and final values of the controllable parameters, and 
not on the path. Mathematically, the thermodynamic 
length is defined by a Riemannian metric on the man- 

Among other 
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ifold of equilibrium ensembles 
useful physical properties, the thermodynamic length 
bounds the dissipation of slow, but finite time trans- 
formations ll|, [l3|. Moreover, the ability to measure 



thermodynamic length and free energy change from out- 
of-equilibrium measurements indicates that these equilib- 
rium properties influence the behavior of driven systems 
even far-from-equilibrium. 

Thermodynamic length was originally defined using 
the second derivatives of a thermodynamic potential with 
respect to its natural variables [tj, [13]. However, this 



definition only works for microscopic systems when the 
controlled variables are intensive (e.g. temperature) [l5| . 
To circumvent this restriction, herein we will redefine 
the thermodynamic length in terms of Fisher informa- 
tion [It], [h|. This approach is equivalent to the orig- 
inal definition for large systems in the thermodynamic 
limit (TBI, [HI], but can also be applied, without restric- 
tion, to microscopic systems, or to problems outside of 
thermodynamics entirely. 

Given a family of probability distributions 7r(a;|A) for 
outcomes x that vary smoothly with a collection of pa- 
rameters A = {A 4 }, the Fisher information matrix [2fj|.[2lT] 
is 



Xjj(A) = / dx n(x\X) 



d\mr(x\A) 91n7r(x|A) 



dX l 



(1) 



The length of a path A(s) for s 6 [0, 1] in parameter 
space measured using the Fisher metric (also known as 
the Fisher-Rao, Rao or entropy differential metric) is [13] 



£ = 



E 



dX^s) 
ds 



2y(A( S )) 



d\'{s) 
ds 



1/2 



ds. 



(2) 



The Fisher matrix 2^ acts as a metric tensor and equips 
the manifold of parameters with a Riemannian met- 
It is also useful to define a related quantity, 



nc 
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the Fisher divergence 



J 



E 



dX^s) 
ds 



2y(A(s)) 



dX j (s) 
ds 



ds 



(3) 



The length and divergence are connected by the relation 
J > C 2 due to the Cauchy-Schwarz inequality. 

The Fisher metric can be applied to any family of prob- 
ability distributions. Here, we focus on probability dis- 
tributions of a system in thermal equilibrium. In the 



2 



canonical ensemble [22|, |23j , the probability of a micro- 
state x is 



7r(ar|A) = exp |/3-F(A) - @E{x, A)} , 



(4) 



where /? = is the inverse temperature T of the 

environment in natural units (fee is the Boltzmann con- 
stant), E(x, A) is the energy of the system, which depends 
both on the internal state x and the external control pa- 
rameters A, and F(X) is the free energy 



0F(\) = -ln^exp{-/3£:(x,A)} . 



(5) 



For a case in which A is a single controllable parameter, 
the Fisher information is 



/(A) 




dE(x,X )' 



(6) 



where (• • • )a indicates an ensemble average over the dis- 
tribution 7r(x|A). 

Let us consider two examples. First, suppose the sys- 
tem under examination is a single polymer, and the pa- 
rameter under control is the end-to-end distance L (Con- 
cretely, an RNA hairpin, with DNA handles, attached to 
beads held by a translating optical trap [3J]). The in- 
stantaneous tension T = ^ is the force exerted on the 
polymer by the apparatus constraining the distance be- 
tween the polymer ends. The Fisher information for this 
system is equal to the variance of the tension at equilib- 



X(L)=/3 2 (((T)-T) 2 ) , 



(7) 



and therefore the thermodynamic length [Eq. |(3J)] is equal 
to the cumulative root-mean-square equilibrium fluctua- 
tions in tension on the molecule. 

On the other hand, suppose control is exerted by ap- 
plying constant tension to the ends of the polymer. The 
total energy is then a linear function of length and ten- 
sion, E(x,T) = U(x) — TL(x), and the Fisher informa- 
tion is equal to the variation of the end-to-end polymer 
length at equilibrium. 



X(T)=/3 2 (((L> 



L) 



(8) 



Again, the Fisher information has a simple physical in- 
terpretation in terms of equilibrium fluctuations, and the 
thermodynamic length [Eq. is equal to the cumula- 
tive root-mean-square fluctuations along the path. If, as 
in the second case, the energy is a linear function of the 
control parameter, the Fisher information is equal to the 
second derivative of the free entropy 15, H, 3], but this 
is not true in the general. 

To demonstrate how to measure thermodynamic 
length in far-from-experiments, we will model the dy- 
namics of the system as a driven, discrete time, inho- 
mogeneous Markov process SIH Hi. The microscopic 



state of the system is denoted by x, and the history of 
the system will be denoted by x a .b = {x a , x a +i, ■ ■ ■ , x b }, 
where a < b. The time reversed trajectory is denoted by 
ib.a = {xb, • • • ,x a }. The internal energy E(x,X) 
depends on a control parameter A that varies according 
to a predetermined protocol, A ab = {A a , A a+ i, • • • , A&}. 
The protocol of the conjugate time reversed experiment 
is Ah, a = {At, A&-i, • • • , A a }. At each integer time t, the 
control value changes from At to Aj+i in the forward pro- 
tocol and At+i to At in the reversed protocol. Between 
these time points A is constant. 

The probability of observing a particular trajectory 
x a .b as the system is driven from thermal equilibrium 
by the protocol K a b can be written as 



b-l 



P a ,b[ X a..b] = ^(x a \X a ) \\Pt(x t +l,X t ) 



(9) 



where 7r(x a |A a ) is the initial equilibrium probability dis- 
tribution with fixed A [Eq. ([I])], and pt(xt+\,Xt) is the 
probability of transitioning from state xt at time t to 
state xt+\ at t + 1, given the value of A at time t. These 
transition probabilities satisfy the stochastic property 



^2pt(x',x) 



(10) 



which ensures conservation of probability, and the bal- 
ance condition 



Pt(x',x)iv(x\\ t ) = 7r(x'|A t ) 



(11) 



which ensures that the canonical ensemble is the station- 
ary distribution of the dynamics. 

The transition probabilities of the forward and time 
reversed dynamics are related, since, at equilibrium, the 
transition x — > x' in the forward dynamics has the same 
probability as the transition x' — * x in the reversed dy- 



namics 



24 l26l ] . given a fixed A. Explicitly, the reverse 



time transitions are related to the forward time transi- 
tions by 



Pt(x,x') ir(x'\\ t ) = Pt(x',x) 7r(x|A t ) . 



(12) 



A direct consequence of this time reversal symmetry is 
the work fluctuation theorem 0, 0, : the ratio of the 
probabilities of the forward and reverse trajectories is the 
exponential of the observed dissipation along the forward 
trajectory. 



Pg,b[ x a,b] 
PbA X b,a] 



f3W a , b lx a , b ]-f3AF atb 



,Va.,b[ x °.,b] 



(13) 



Here AF a ,b = Fb — F a is the free energy change and 
6-1 

W a , b [x a!b ] = [E(x t+1 ,Xt+i) - E(x t+1 , A t )] (14) 

t — a 
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is the work transferee! to the system during the for- 
ward process [2f| |27|. The dissipation T> a ^[x a ,b] = 
(i(Wa,b[ x a,b\ ~ AF^), the irreversible increase in en- 
tropy along the forward trajectory, is proportional to 
the difference between the work and the free energy 
change. Note that work, free energy change and dissi- 
pation are all odd functionals under a time reversal, e.g. 

Wb, a [%b,a\ = -W a ,b[x a ,b]- 

We can express the trajectory ensemble average of an 
arbitrary trajectory dependent function ^[xa^b], starting 
from thermal equilibrium, as 



(F[ x a,b])ab = E P aA x a,b] F[x a ,b\ 



(15) 



and similarly for the conjugate process 

(F[x bt a]) b a = ^ Pb,a[xb,a] F[Xb,a] (16) 

in which !F[Xb a ] = !F[x a .b] is an even functional under 
time reversal. 

A key result in our development links two different 
trajectory ensemble averages 



F[x a ,b] 



a . b 



.-T>o 



F[x a ,b] 



0,T 



(17) 



where < a < b < T. Given a protocol Ao,t, we can 
extract the value of a trajectory ensemble average over a 
subinterval h- a ,b, as if the system began in equilibrium at 
an intermediate time a, by re- weighting the observations 
by the exponential of the dissipation from the initial to 
intermediate time. 

This result follows directly from the work fluctuation 
relation [Eq. (Tl3|) ]. and the Markovian property of the 
dynamics. 



e -V . a [*o. a \ T[x b] ) = / e -V , aMj r [x A (18) 
0,T \ / 0,6 



-T>h, a ]x b 



^[Xb,a_ 



6.0 



{HXa.b]) a , 



We truncate the time interval of the trajectory ensemble 
average using the stochastic property in Eq. (fTO]) . ap- 
ply a time reversal with the work fluctuation theorem in 
Eq. (fT3|) . truncate again, and apply a second time rever- 
sal. This result generalizes previous trajectory ensemble 
averages of Hummer and Szabo [3] and Chelli et al. [28| . 

We can now use this relation to extract the thermo- 
dynamic length from far-from-equilibrium experiments. 
The discrete time analogs of the Fisher length and diver- 
gence are the cumulative Jensen-Shannon length 



is 



V8 



T-l 

E 

t=0 



JS 



(7r(s|At),7r(a:|At+i)) (19) 



and cumulative Jensen-Shannon divergence [if 



T-l 



J JS = 8 JS(7r(i|A t ),7r(i|A t+ i)) 



(20) 



t=o 



Here, JS(j?i, p%) is the Jensen-Shannon divergence [29ll3 
between two probability distributions p\ and p 2 , 



JS(pi,P2) = -^pi(x)ln T 

X 2 

+ \^P2{x) In j 



Piix) 



Pl(x) +P2(x)} 
P2(x) 



(21) 



[pi(x) +P2{x)] 



The Jensen-Shannon length is less than the Fisher length 
£js < A and approaches equality as the step size along 
the path decreases (l5| . 

We can use the contracted trajectory average Eq. (fTTjl . 
and the canonical probabilities [Eq. @\ to write the 
the Jensen-Shannon divergence between any pair of time 
points along the path 



(22) 



o,t+i 



1 + e- v t+i-.t 



TA 



as a trajectory average of the dissipation T> along the for- 
ward and reverse protocols. While JS(-7r t , TTt+i) is defined 
in terms of averages over equilibrium probability distri- 
butions, it can be related to trajectory ensemble averages 
of processes driven arbitrary far-from-equilibrium. The 
derivation of Eq. (|22[) requires the time reversal symme- 
try in Eq. (fHf . 

We now encounter an apparent complication. The dis- 
sipation T> a £ — (3 {W a ,b — ^F a ,b) depends on both the 
work and the free energy. Therefore, we must also de- 
termine the potential of mean force, the free energy as 
a function of A, along the entire path. This problem 
of extracting free energy profiles from out-of-equilibrium 
work measurements (rather than just the difference in 
free energy between the initial and final ensembles) has 
attracted recent attention [28L l3lL l32t| . Here, we will solve 
this problem by adapting Bennett's maximum likelihood 
method 0, [33|, [34], [35|, which, as we shall see, is inti- 
mately linked to the thermodynamic divergence [l5T | . 

Suppose we have taken measurements of the work dur- 
ing N repetitions of a protocol A a ,& and another N mea- 
surements from the conjugate protocol A(, )0 . Each rep- 
etition begins in thermal equilibrium with the control 
parameter fixed at A Q or A&. Then the Bennett log- 
likelihood that the free energy change AF a ^ has a par- 
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ticular value is 15 1 



£(AF a , b ) =J] In 



1 



n=l l + e 
N 



1 



(23) 



+ Vln 

~l 1 + e-^'+pAF^ 

where Wy 1 } and are the work measured between a 

and during the nth repetition of the forward and reverse 
experiment respectively. 

Next we extend this result using the contracted trajec- 
tory average [Eq. (|17[) ] to estimate the likelihood of the 
free energy change between any two points along the pro- 
tocol. In particular, we can estimate the log-likelihood 
for the entire free energy profile by summing the log- 
likelihood for every pair of neighboring time points 



T-l 



*m) - E 



t=0 



N 



71 = 1 



]T e-W In 
l 

N 

E 



i 



1 + e^M+i 
1 



(n) 



(24) 



l + e *+!■' 



Here, {F} = (Fq, F\, ■ ■ ■ , Ft) is the free energy profile. 
Since only differences in free energy are relevant, one 
free energy is set to zero or some other convenient ref- 
erence. Since each experimental realization of the for- 
ward and reverse protocols begins in equilibrium, the 
probability of each forward and reverse realization is 
Pq,n[xo,n] and Pn,o[xn,o] respectively. Hence, we have 
written this expression using the measured dissipation 
V^l = 0(wj$ - AF Qib ) so that the relationship with 
Eq. (|2"2"|) is clear. To within an additive constant, the total 
Bennett log-likelihood is proportional to the cumulative 
Jensen-Shannon divergence. Therefore, we can simulta- 
neously determine the potential of mean force, the ther- 
modynamic divergence, and the thermodynamic length 
from the same collection of work measurements. We first 
determine the free energy profile {F} that maximizes the 
log-likelihood, which immediately provides an estimate of 
the thermodynamic divergence, 



TN\n2 



(25) 



We can then calculate a maximum likelihood estimate of 
the thermodynamics length £jg in Eq. (|19p using the free 
energy profile. 

In principle, it should be possible to measure thermo- 
dynamic length in single molecule experiments. One cap- 
tures a RNA hairpin in an optical tweezer, and repeat- 
edly measure the force on the RNA molecule as a function 
of extension. These far-from-cquilibrium force extension 
curves will then yield the potential of mean force, plus 
the thermodynamic length and divergence of the proto- 
col. Far-from-equilibrium measurements of thermody- 
namic length have interesting implications for nano-scale 



machines and biological motors, since the square of this 
length bounds the dissipation during finite time proto- 



cols llj. Hence, thermodynamic length is intimately 



connected with the useful work that a system can per- 
form. It would, for example, be interesting to measure 
the thermodynamic length along the cycle of a molecular 
motor. For a machine operating at a finite rate, but oth- 
erwise optimized to minimize dissipation, it is expected 
that the rate of change in thermodynamic length along 
the cycle would be constant flij ]. 
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